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We present a general Lagrangian formalism that allows the treatment of vorticity. We 
give solutions for the rotational perturbations up to the third-order in a flat background 
universe. We show how the primordial vorticity affects the evolution of the density fluctuation 
in high-density regions. 

§1. Introduction 

The study of the galaxy formation based on gravitational instability has not 
revealed the origin of the vorticity of galaxies. The possibility that galactic spins 
are the relics of primordial vorticity has not been taken into account. Because the 
standard linear Eulerian theory (see, for example, Peebles' bookliP for the detail 
of the Eulerian theory) shows that the vorticity decays at late time and does not 
couple to the density enhancement. For this reason, the methods that generate the 
vorticity posteriorly are investigated. DoroshkevichtP explained the acquisition of 
the vorticity by the formation of ^hock fronts in the pancake. The explanation by 
tidal torques was made by Hoyle.™ 

However, can we actually ignore the effects of the primordial vorticity in the 
galaxy formation process? To answer this question, we must reconsider the validity 
of the Eulerian theory. In the Eulerian treatment, density fluctuation is perturbative 
quantity, so it must be small. Therefore, the Eulerian approach is limited to the 
description of weak density fluctuation. For this reason, to discuss the behavior of 
the vorticity in the high-density regions, we must prepare the theory beyond the 
Eulerian theory. The Lagrangian theory improves such shortcomings of the Eulerian 
theory. The most advantageous point of the Lagrangian theory is that the density 
fluctuation is non-perturbative quantity. Therefore, it is unnecessary that the density 
fluctuation is expanded in a perturbation series and that we impose the condition 
of the smallness of the density fluctuation. By using the property of the Lagrangian 
theory, we can investigate the behavior of the primordial vorticity. BuchertQ^ solved 
the first-order Lagrangian perturbative equations and showed that the primordial 
vorticity was amplified in proportion to the enhancement of the density fluctuation 
by deriving the relation between the vorticity and the density fluctuation. Barrow 
and SaichcP solved the vorticity equation on the inhomogeneous background and 
pointed out that the growing directions of the vorticity were lying in the plane of the 
pancake. These analyses based on the Lagrangian theory shows that the primordial 
vorticity is not negligible in the high-density regions such as the pancake. 

In this paper, for the purpose of preparing a tool to treat the behavior of the 
vorticity in the high-density regions, we extend the previous works and give solutions 
including the vorticity for the Lagrangian perturbative equations up to the third- 
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order. Furthermore, we obtain the relation between the vorticity and the density 
fluctuation. The reason we proceed to the calculation up to the third-order is that 
the third-order solutions will cover the main effects since the density fluctuation is 
described by the determinant of a 3 x 3 matrix. In the Lagrangian approach, we 
assume that background is Friedmann-Lemaitre-Robertson- Walker (FLRW) universe 
with no cosmological constant and curvature. There is no definitive observational 
evidence that our universe is actually flat. We apply this assumption for simplifying 
calculations. 

This paper is organized as follows. In §|2], we summarize the basic formulae of the 
Lagrangian theory. In §||, we present the perturbative approach to the Lagrangian 
equations including the treatment of the vorticity. In §||, we summarize the main 
results of this paper. 

§2. Basic equations in Lagrangian form 

In this section, we give a summary of the Lagrangian theory. We particularly 
focus on the treatment of the vorticity in the Lagrangian theory. The detail of the 
Lagrangian theory is found, for example, in Ref. ||). 

2.1. Density fluctuation in Lagrangian form 

In this subsection, we focus on the expression of the density fluctuation in La- 
grangian form. 

In the Eulerian perturbation theory, the density fluctuation 5 is a perturbative 
quantity. On the other hand, in the Lagrangian form, the density fluctuation is given 
in the formally exact form as follows: 

*-J-'-l,J=-d-(g, (2-1) 

where x are expanding coordinates in the Eulerian space and q are the Lagrangian 
coordinates which are defined by the initial Eulerian position, namely, q := x(t = to). 
(The physical distance is obtained according to the law r = a(t)x, a(t) being the 
scale factor.) The determinant J of the Jacobian of the transformation from x 
to q is basic quantity to represent the density fluctuation. Since J is defined by 
the Lagrangian partial derivative, we must evaluate x in the Lagrangian space: 
x = x(q,t). Therefore, it is necessary to clarify the relation between the Eulerian 
coordinates and the Lagrangian ones. The displacement vector is the quantity that 
relates x with q. This quantity represents time evolution from the initial position q. 
By using the displacement vector A(q, t), we can write the relation between x and 
q as x = q + A(q, t). If we obtain the concrete form of A, the density fluctuation is 
calculated from Eq. fl2T|) . Since Eq. ( |2-1|) is derived from mass conservation between 
the Eulerian and the Lagrangian space, this equation is hold irrespective of the 
existence of the vorticity. 

Now we shall consider how the effects of the vorticity on the density fluctuation 
are reflected. A conventional way to investigate the effects of the vorticity is that 
we split A into the longitudinal mode and the transverse one with respect to the 
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the Eulerian coordinates and then focus on the behavior of the transverse mode. 
However, since A is evaluated in the Lagrangian coordinates, we shall perform this 
decomposition with respect to the Lagrangian coordinates: 



AW + A\ V x A 11 = 0, V» • A J 



0. 



(2-2) 



where the superscript || and _L denote the longitudinal and the transverse mode with 
respect to the Lagrangian coordinates, respectively. Although the quantity A^ does 
not correspond with the vorticity, we shall express the density fluctuation by using 
A" and A^. The method to express the vorticity in the Lagrangian space will be 



presented in next subsection. From the definition of J in Eq. (2-1), we obtain 
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and the comma denotes the partial derivative with respect to the Lagrangian coordi- 
nates. By considering the transverse mode A" 1- , we can investigate how the vorticity 
contributes to the density enhancement. 

In next subsection, we derive the Lagrangian equations for A" and A ± which 
are necessary for the calculation of the density fluctuation. 

2.2. Equations for longitudinal and transverse mode in Lagrangian coordinates 

As mentioned in previous subsection, we derive the Lagrangian equations for 
A" and A 1 - here. ^The derivation of the equations is based on the methods by 
Buchert and G6tz,LP and KasaitP. Furthermore, we also present the expression of 
the vorticity in the Lagrangian space. 

First, we shall derive the Lagrangian equations. We denote the density of dust, 
the background density, the peculiar velocity and the peculiar gravitational acceler- 
ation by p, pb, v and g, respectively. The peculiar velocity given by the form v = ax 
represents the deviation from the uniform Hubble flow. Since we assume the FLRW 
universe, the scale factor a is proportional to t 2 ^ 3 . We consider that a = i 2//3 . The 
basic equations in Newtonian cosmology that these quantities obey are 



dp 'a 1 , . 
-£ + 3-p+-X7 x -(pv) 
ot a a 



0, 



(2-5a) 



4 



M. Sasaki and M. Kasai 



dv a 1 , _ , ,„ „. . 

- 5 r + ~« + -(wV a .)t>=0, 2-5b 
at a a 

V I xg = 0, (2-5c) 
V* • = -47rGa(p - p 6 ), (2-5d) 
where V x is the Eulerian nabla operator. We obtain from Eqs. ( |2-5a| ) and ( j2-5d|) 



v., 

and from Eqs. ( 2-5b| ) and ( |2 5<j ) 

V, X 



-£ + 2-g- 4nGpv 
ot a 



0. 



dv a 1 . _ . 
— + -v + - (v ■ V x )v 
at a a 



0. 



(2-6a) 



(2-6b) 



Eqs. ( |2-6a ) and ( |2-6b ) are the equations for the longitudinal and the transverse mode 
with respect to the Eulerian coordinates, respectively. We rewrite the system ( |2-6| ) 
to the Lagrangian equations by transforming the Eulerian partial derivative into the 
Lagrangian one 
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and introducing the Lagrangian time derivative 
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(2-9a) 

%jj ~ -eijkBijV{A k ,i), (2-9b) 
where T> is the differential operator satisfying 



V(e ijk Ai 



m/ri , (d 2 F adF\ 



(2.10) 



Cijk is Levi-Civita asymmetric tensor, ei23 := 1, and V g is the Lagrangian nabla 
operator. The detail of the second term in right-hand side of (2-9a) is in Appendix 
[A|. Eqs. (2-9a) and ( 2-9b[) are the equations for the longitudinal and transverse mode 
with respect to the Lagrangian coordinates, respectively. The system fl2-9|) is the 
Lagrangian equations which are necessary for obtaining the density fluctuation. 
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Next, we shall consider the expression of the vorticity in the Lagrangian space. 
Here, we define the vorticity as w := (l/a)V x . x v. By transforming x into q, we 
obtain 

1.1 , 
Wi = - V„ x v \i + -eijkBijVkj. (2-11) 
a a 

However, there exists another expression of the vorticity. In the Lagrangian repre- 
sentation, the equation for the transverse mode with respect to the Eulerian space 
( 2-6b| ) is solved exactly as (see Ref. f|)) 



w = -1(1 + 6) (w(q, t ) • V,) x(q, t), (2-12) 

where w(q,t$) is the primordial vorticity. Eq. (|2-12j ) represents that the vorticity 
couples to the enhancement of the density fluctuation. This fact is not respected 
in the Eulerian theory. Naturally, these two expressions of the vorticity ( 2-lip and 



( 2-12| ) are equivalent. We stress that the Lagrangian method enables not only to 



calculate the density fluctuation straightforward but to clarify vorticity-density re- 
lation. 

In next section, we will solve the Lagrangian equations by perturbative method 
and obtain the concrete expression of the density fluctuation and the vorticity. 

§3. Lagrangian perturbative approximation 



In this section, we solve the Lagrangian equations (2-9a) and ( |2-9b ) by perturba- 



tive method up to the third-order and obtain the determinant J. By calculating J, 
we obtain the concrete expression of the density fluctuation and the vorticity from 
Eqs. ( |2-1| ) and (|2T2 ), respectively. Our treatment is different from the standard 



Lagrangian perturbation theory in considering the vorticity. Naturally, when the 
irrotational condition is assumed, our results are consistent with the standard La- 
grangian perturbation theory. The reason the perturbative approach is introduced 
is that the right-hand sides of the system (|2-S| ) are non-linear. However, since the 
density fluctuation itself has no restriction, Lagrangian perturbative approximation 
enables to describe the high-density regions. (The restriction by the Lagrangian per- 
turbation approach is that we require the smallness of the deviation of the particle 
trajectory from the uniform Hubble flow.) In the Lagrangian perturbation approach, 
the perturbative quantity is the displacement vector: 

x(q,t) = q + A (1) (q, t) + A (2) (q, t) + A (3) (q, t) + • • • , (3-1) 

where Am (q, t) corresponds to the first-order approximation, A( 2 ) (q, t) to the second- 
order approximation, and so on. We derive the perturbative equations for A^), 
A( 2 ) and A( 3 ) and then solve these equations. Hereafter, we call the longitudi- 
nal (transverse) mode with respect to the Lagrangian coordinates simply "the lon- 
gitudinal (transverse) mode". 

3.1. First-order solutions 

The first-order solutions for the longitudinal and the transverse mode are ob- 
tained by Buchert.& 
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First, we consider the longitudinal mode. The equation for the first-order longi- 
tudinal mode is 



d_ 

df 



a 3 {p(A[ l 1) )-47rG /0f) A[ l 1) } 



V 



(3-2) 



where K is an arbitrary vector potential. The existence of the quantity K means 
that the decomposition such as A = A' + A 1 - is not unique; we may always add to 
A" a rotation-free solution of V • A" =0. We can impose additional condition for 
eliminating this freedom. So we shall assume V ? 
obtain 

H C dC 



K = 0. Integrating Eq. (3^2), we 



Mi)) 
By choosing C = 0, we obtain 

A 11 - 



AirGpbAfa = — , — = 0. 



dt 



(3-3) 



(3-4) 



(The Lagrangian perturbative solutions at any order separate with respect to q and 
t. See Ref. ||).) The terms proportional to i 2//3 and t~ 1 are called the growing mode 
and the decaying mode, respectively. 

Next, we consider the transverse mode. The equation for the first-order trans- 
verse mode is 

V(Afq) = V q <p, (3-5) 

where ip is an arbitrary scalar potential. We assume = for eliminating the 
freedom. The solution is 

Afc^t-Wxfaq). (3-6) 

We shall call this solution the rotational mode. 

3.2. Second-order solutions 

The general irrotational second-order solutions are obtained by Buchert and 
Ehlers& in the case of density parameter [2 = 1. The leading terms of the irrota- 
tional second-order solution in the case of arbitrary density parameter are derived by 
Bouchet et al..H2P In order to avoid notational complexity, we neglect the decaying 
mode -X"ll « hereafter. 



First, we consider the longitudinal mode. The equation for the second-order 
longitudinal mode is 

|[a 3 {p(V 9 .A[ l 2) )-4^Gp 6 V f/ .Al | 2) } 



-B 



(i)v 



The particular solution is 



^ 2) =t^l 2a) (q)+t^xl b) (q), 



"(2a) 



(3-7) 



(3-8) 
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where 

, 2 



V 9 "A (2a) " 14 



(1) 



(3-9) 



V 9 ' X \2b) ~ \ X {l)id X t)i,i- ( 3 " 10 ) 

Although the solution for Eq. (3-7) consists of the homogeneous and the particular 
solution, we omit the homogeneous solution. We focus on the behavior of the first- 
order quantities, namely, the growing mode x\^ and the rotational mode X^y m 
higher-order. Since the homogeneous equation is the same form as the first-order 
equation, the behavior of the homogeneous solutions is investigated easily. Now the 
solution ( |3-8| ) consists of two parts. One is -^[' 2a ) described by using two growing 

modes. Another is ■X'Lm described by using two rotational modes. In the case where 

the rotational mode is considered, the solution AIL includes the Xj~^ , although AIL 
describes the longitudinal mode. Furthermore, there is no cross term consisting of 
an d Xj^y This cross term is included the transverse mode as mentioned below. 
Next, we consider the transverse mode. The equation for the second-order trans- 
verse mode is 

rD (^ijk^f 2 )k,j) = - ( -ijkB{i)i j V(A { i) kl ). (3-11) 

In the same sense as that of the longitudinal mode, we shall obtain the particular 
solution for Eq. ( |3T1| ). The particular solution is 

Af 2) =t l /*Xf 2y (3-12) 

where 

% * ^(2) 



: -3eijkX^jjX( X y )k . (3-13) 

The second-order transverse mode is induced by the coupling of the rotational mode 
to the growing mode. Note that this coupling term is growing. On the other hand, 
the cross term of two -X"Pn or two X^ does not exist in Eq. (3-12). Therefore, if 
the rotational mode is zero, the second-order transverse mode vanishes. Actually, we 
must consider the fact that the decaying mode couples to X^L, or X^y However, 
the coupling of the decaying mode is different from the coupling of X/^ to the 

in decaying at late time. Although both the decaying mode and the rotational 
mode decay, the former is more effective when they couple to the growing mode in 
higher-order. 

3.3. Third-order solutions 

The irrotational third-ordeii solution is obtained, for example, by Buchert0', 
Bouchet et al.,ll3> and CatelanE^P. 

First, we consider the longitudinal mode. The equation for the third-order lon- 
gitudinal mode is 



d_ 

dt 



[a 3 {p(V ff -A^ ) )-47rG P6 V,-^ ) } 



8 



M. Sasaki and M. Kasai 



R d 



a z [v{A {l)j)i ) -AnG Pb A mi ] 
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The particular solution is 
A 

where 



l (3) 
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(3-17) 
(3-18) 



The third-order solution (3-15) does not contain the term consisting of two X^ 



and one -^m- Furthermore, note that the third-order solution includes the term 

proportional to i _1 lni, although the solutions up to the second-order are described 
in the form of simple power of t. Since this term consists of three X^, the term is 

anticipated proportional to namely, (t -1 / 3 ) 3 . However, the solution proportional 
to t^ 1 is same form as the first-order decaying solution, that is, the homogeneous 
solution of the Lagrangian equation. Therefore, the term proportional to hit is 
necessary. 

Next, we consider the transverse mode. The equation for the third-order trans- 
verse mode is 



^( e ijk A (3)k,j) — ~ e ijk 
The particular solution is 
X 

where 
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(3-20) 

(3-21) 
(3-22) 
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X (l)l,j X (2)l,k 
1 



2X\L, -X. 



(l)l,j^(2b)k,l 



(3-23) 



£ yA;^(3d)fc,j - 3 e yfc^(t)/j^(2fe)fc,/- (3-24) 

Note that even though the irrotational condition is assumed, the non-zero transverse 
mode exists in the first term in the right-hand side of Eq. ( |3-20| ), as BuchertE^ 
pointed out. In deriving ( |3-21| )~(3-24), by using the following relation 

H r„ii „i „n „i i 1 



X (l)rn,l X {l)l,k X (l)k,l X {l)l,m 



3^*^(1)^,1 



^(2)J,fc ~~ ^(2)fc,« 



€ ijk X (i)j t l X (l)l,m X (l)m,k 



3 e yfe^(i)i,j 



-^(2)Z,fc ~~ -^(2)fe,i 



(3-25) 
(3-26) 



we write the triplet of -X"(i) as the term including the second-order quantity. By the 
same reason in the case of the longitudinal solution, the term proportional to hit 
exists in the third-order transverse solution. 



3.4. Relation between vorticity and density fluctuation 

In this subsection, we calculate the concrete form of the density fluctuation 
and the vorticity by using the Lagrangian perturbation formalism we have obtained 
above. 

From the definition of the density fluctuation Eq. ( |2-1| ), we obtain the density 
fluctuation as follows: 

1 



det 



kj + + A (2 )ij + A {3)iJ 

1 



1 



1 + J(i) + J( 2 ) + J(3) + 



where 



Jm = i 2/3 V„ • X 



(i)' 



(3-27) 
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Jt 
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: t 4/ 3 



(y q -x 



x (i)ij x (i)j,i 



/1/3V-II v± , ^-2/3yl Y ± 
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(3) 



q { X l)i,i X la) hj - X (i)i,j X (L)j,J + l det ( X . 



+ t 
+ 



3 X (l)i,i X (2)j,i + ^ X (2a)i,j yY (l)j,i 



■XT 



-X(i)ijX {2)jA + 



{ij'.l 
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II 



+ _r i (201nt-63)^ 6) . i X ( - )j . i 
+ ^t- 1 (19-101ni)X ( -L ) . i X ( i ) . fe X 



(i)M' 



(3-30) 
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Although we represent J only up to the third-order here, we can calculate J up to the 
ninth-order. Recall that it is unnecessary that the density fluctuation is expanded 
perturbatively in spite of J being the perturbative quantity. By linearizing Eq. (3-27), 
we recover the Eulerian linear behavior presenting that the vorticity decouples from 
the density enhancement. From the quantities J( 2 ) and J( 3 ), we see that the coupling 

of the rotational mode X^ to the growing mode -XjL plays an important role in 
the enhancement of the density fluctuation. Because this coupling is growing. The 
decaying mode we have omitted also can contribute to the enhancement of the density 
fluctuation. If we write the term containing one X^L * and two xjL in J^, this 

term is proportional to t 1 / 3 , namely, i _1 (i 2 / 3 ) 2 . However, since the term coupling of 
the rotational mode is proportional to t, the rotational mode makes more important 
contribution to the enhancement of the density fluctuation than the decaying mode. 
The importance of the coupling of the rotational mode to the growing mode is also 
seen from vorticity-density relation (see Eq. ( |2T2D ) 

1 (w(q, t ) ■ V q ) (q + A (1) + A (2) + A (3) ) 



w 



1 + J(l) + J(2) + J(7t\ + ■ ■ ■ 



(3) 



(3-31) 



The vorticity-density relation shows that the behavior of the vorticity depends on 
both the density fluctuation (<5 or J) and the deformation of the fluid (x{j or A^j). 
For realizing the concrete form of the vorticity, we shall expand Eq. (|3-3l| ) pertur- 
batively. (Note that the vorticity is actually the non-perturbative quantity in the 
Lagrangian form.) We obtain 



w 
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(3-32) 
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+t -2/3 |-I W(1) (g)X ( -L )aj6 X ( i )6ia + • V,) X[ 

Here, we assume 



(26) 



(3-34) 



w(q,t ) := tt>(i)(g) 



--V D x X 



(!)• 



(3-35) 



Same results are also derived from Eq. (|2-ll| ), that is, the transformation of the 
Eulerian coordinates to the Lagrangian one (see Appendix [C]) . The results ( |3-32| ) 
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is the same as the behavior of the vorticity in Eulerian description. Furthermore, 
from the results (3-33) and (3-34), we see that the amplitude of induced vorticity 
becomes strong by the coupling of the rotational mode (initial vorticity) to many 
growing modes. This shows that the primordial vorticity is effective in the high- 
density regions. The exact calculation of the vorticity Eq. (3-31) will confirm this 
fact. 



§4. Summary and conclusion 



In this paper, we have investigated the relation between the vorticity and the 
density fluctuation by preparing the rotational Lagrangian perturbation theory. The 
fact that the vorticity couples to the density enhancement is already discussed by 
Buchert.u) However, we obtain more concrete vorticity-density relation by solving 
the perturbative equations up to the third-order for both longitudinal and transverse 
displacement vector. The coupling of the transverse mode to the growing mode 
enables the vorticity to contribute to the enhancement of the density fluctuation. 
Furthermore, growing vortical modes are induced by the coupling of initial vorticity 
to the growing mode. These results show the fact that the primordial vorticity is not 
negligible in high-density regions. Therefore, we must reconsider the importance of 
the primordial vorticity effect on the density fluctuation. Now we also have pointed 
out the effects of the decaying mode. It is possible that the first-order decaying 
mode contributes to the enhancement of the density fluctuation by coupling to the 
growing mode. However, since the transverse mode is proportional to i -1 / 3 and the 
decaying mode to t , the former makes more important role in the enhancement of 
the density fluctuation. 

The results we obtained in this paper give a powerful tool to treat the high- 
density regions and behavior of the vorticity. It is of interest to investigate quanti- 
tatively how the vorticity affects the non-linear evolution of the density fluctuations 
in the universe. It will be the subject of future investigation. 
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Appendix A 

Detail of Lagrangian equation for longitudinal mode 

For the perturbative calculation, we give the detail of the second term of the 
right-hand side of the longitudinal equation (2-9a) as follows: 



V x ■ [[A ■ V x ) V(A) - (V* • V{A)) A 



12 
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+B ij A jtk V(A k , i ) + B tj V(A j:k )A kA - BijAj^Ay) - B^A^A^ 
+Bij :k AjV(A kji ) - Bij tk A k T>(Aj t i) 



+Bi 



B kl A t V{A jyk ) - B kl V(Ai ik )Aj 



Appendix B 

Calculation of B{j 



(Al) 



When we solve the perturbative equations, the quantities and B( 2 )ij are 



necessary. So we present the derivation of their quantities here. 

The quantity B^ is defined by using the cofactor C{j as follows: 



B 



1.1 



71=1 

Using Eq. flB-l| ), we obtain 



(») 



£jab£icd%a,c%b,d' 



B, 



B 



(2)V 



^(l)i,fc^(l)fc,j -^(2)iJ- 



(B-l) 

(B-2) 
(B-3) 



Appendix C 

Useful relations to obtain the expression of the vorticity 



We have presented the expression of the vorticity in Eqs. ( |3-32D, (3-33) and 
(3-34). These results are derived by vorticity-density relation Eq. ( |2T2| ). We also 
obtain the same results from Eq. ( 2-llj ), that is, the transformation of the Eule- 
rian coordinates to the Lagrangian one. We list the useful relations to obtain the 
expression of the vorticity from Eq. ( 2T1 ) as follows: 



e lmnX(i)n tm .X (lfti — € ijk X (i)j t i X (i)l t k> 



,x 



(l)n,rn yl (l)i t l ~ c «7 



^ijk 



X. 



ijk^(l)k,j X (l)l,m X (l)m,l 



(1) 



X? 



x 



(i)k,i 



einhX, 



_L 



X(-,\A lX, 



ijk^(l)j,l^(l)l,m X (l)m,ki 



e irrm^(i) rai m^(2)i,I — e ijk 



x 7 . 



x, 



X 



X, 



(C-l) 

(C-2) 
(C-3) 
(C4) 



e ijk X (l)l,j 



X 

II 



(2)/,fc 



(2)fc,Z 



3e ^fc X fl)/,m X (l)i,7 X (l 



3eiofcXjL • ,X. 



X? 



{ i k (l)l,m A (l)m,lfe 



3e 



^(i)j,m^(l~)/,m^(l>,/ 



(l)fc,i^-(l)/,m 
3ej mn -^(i)n,m-^'(l)i,J ' ^-(1)) ' 



^(l)m,i + X (l)l, m X {l)m,j X {i)k > l 



(C-5) 



Non-linear vorticity- density coupling in Lagrangian dynamics 
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^ijkX^ 



x, 



_L 



(2)1, k 



x, 



€lmnX l 



X 



(l)k,l 
x' 



II 

(l)l,m~^(l)m,k 



X. 



Xt 



(l)k,m (l)m,l 



X. 



mn vi -(l)n,m^(l)i,/ 



V 9 • X {1) 



+ tijkX 



,x 



II 



X 



l ' jk (I) 3,1 v (l)l,m yv (l)k,m 



jkX (1 



( V <? " X (l)) + X (i)i, m X (l)m,i 



(C-6) 
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